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AfiSTRACT : 

An  automated  method  for  finding  and  extracting  water  regions  from  Synthetic  Aperture  Radar  (SAR)  imagery 
is  presented.  Three  aucomated  water  finders  are  implemented  for  experimentation  with  two  secs  of  SAR 
imagery  obtained  from  the  areas  of  Elizabeth  City  and  Asheboro,  North  Carolina.  The  input  radar  image 
Is  first  segmented  into  four  categories  of  terrain  features.  They  are  water,  field,  forests,  and  built- 
up  areas.  Water  regions  from  the  segmented  Image  are  extracted  and  the  remaining  image  is  eliminated. 
Noise  appearing  on  both  water  and  non-water  regions  is  smooched  and  thus  completes  Che  automated  water 
finding  cask.  A  second  water  finder  with  a  relatively  higher  threshold  for  the  region  growing  process 
was  designed  for  extracting  shallow  water  regions  from  some  SAR  Images.  Discrimination  of  non-water 
features  having  water  like  gray  values  is  accomplished  by  using  a  third  water  finder  equipped  with  an 
additional  capability  of  computing  and  examining  region  properties. 

KEY  WORDS:  Image  Analysis,  Image  Interpretation,  Image  Processing,  Machine  Vision,  SAR 


INTRODUCTION  SYSTEM  DESCRIPTION 


•  The  problem  with  automatic  extraction  of  water 
bodies  from  radar  imagery  has  been  the  subject 
of  research  for  some  time.  In  the  paec,  various 
pattern  classification  methods  have  been  applied 
to  samples  of  both  Synthetic  Aperture  Radar  (SAR) 
and  aerial  photographic  imagery.  In  many  cases, 
when  the  samples  of  the  water  regions  were  homo¬ 
geneous  and  each  sample  contained  only  the  water 
category,  a  highly  successful  classification 
accuracy  was  obtained  <Fox  and  Chen.  1988).  For 
ell  paec  exporimcnCs  mcncloncd  .tbovc,  an  Image 
sample  size  of  32  by  32  pixels  was  used.  This 
sample  window  size  was  found  Co  be  optimum  for 
these  studies,  even  though  the  size  of  the  sam¬ 
ple  window  can  be  varied  arbitrarily. 

One  of  Che  major  drawbacks  of  a  statistical  pat¬ 
tern  classification  system  is  its  inability  to 
classify  samples  containing  two  or  more  terrain 
categories.  This  problem  occurred  most  often 
when  the  sample  window  was  moved  across  the 
boundary  of  different  image  regions  or  catego¬ 
ries.  A  totally  different  approach  was  proposed 
by  Chen  and  Hevenor  in  a  recent  paper  (Chen  and 
Hevenor,  1990).  In  that  paper,  they  tried  to 
overcome  this  drawback  by  first  segmenting  the 
entire  SAR  image  into  four  categories  of  terrain 
features  where  the  boundaries  were  well  preserv¬ 
ed.  Extraction  or  classification  of  the  terrain 
features  was  then  easily  performed  afterwards. 

In  this  paper,  the  technique  described  In  (Chen 
and  Hevenor »  1990)  was  used  as  a  preprocessor 
to  segment  a  given  SAR  image  into  four  terrain 
categories  of  water,  fields,  forests,  and 
built-up  areas.  The  water  regions  are  then 
extracted,  and  the  remaining  image  is  eliminated. 
The  extracted  water  regions  are  smoothed  and  thus 
completes  the  desired  automated  water  finding 
task.  A  second  water  finder  was  designed  for 
finding  shallow  water  regions  that  requires  a 
relatively  higher  threshold  value  on  the 
region  growing  operation.  A  third  special 
water  finder  was  implemented  for  extracting 
only  the  largest  water  body  from  all  water  reg¬ 
ions  found  from  an  Image.  Discrimination  of 
non-water  features  having  water  like  gray 
values  was  also  accomplished  by  the  use  of  this 
third  water  finder. 


Two  sets  of  SAR  Imagery  were  used  for  testing 
three  automated  water  finders  described  below. 
The  first  sec  of  SAR  Imagery  was  X-baod,  and 
HH  polarization  which  was  taken  over  the  Eliza¬ 
beth  City,  North  Carolina  area  with  a  UPD-4 
radar  system.  Approximately  200  SAR  images 
were  digitized,  and  scored  for  various  feature 
extraction  tasks.  Each  digitized  image  con¬ 
sists  of  312  by  512  plxele  and  each  pixel  con¬ 
tains  8  bits.  Each  image  represents  a  ground 
area  of  approximately  1.6  by  1.6  square  miles. 
The  second  sec  of  SAR  imagery  was  also  X-band, 
but  was  VV  polarization.  This  set  of  imagery 
was  taken  over  Che  Asheboro,  North  Carolina 
area.  Each  Image  also  consists  of  312  by  312 
pixels  of  8-blC  length,  and  an  entire  image 
covers  s  ground  area  of  approximately  0.8 
by  0.8  square  miles.  Two  SAR  images  that  con¬ 
tain  water  regions  were  selected  from  this  set 
of  imagery  for  testing  the  automated  water 
finders . 

Most  of  Che  images  selected  contain  area 
features  such  as  water,  fields,  forests,  built- 
up  areas,  and  some  linear  features  such  as 
bridges,  roads,  railroads,  and  boundaries 
between  area  regions.  The  required  algorithms 
for  performing  automatic  water  finding  task 
were  all  written  In  the  C  programming  language 
and  Implemented  on  a  SUN  4/330  microcomputer 
system  for  the  purpose  of  transferring  them  to 
a  development  laboratory. 

METHODOLOGY 

There  are  several  algorithms  chat  are  required 
CO  automatically  segment  and  extract  water 
regions  from  SAR  imagery.  The  process  for  the 
automated  water  finder  for  SAR  images  consists 
of  the  following  steps: 

1.  Load  the  desired  digitized  SAR  image  from 
the  disk  to  the  computer.  The  image  will 
be  shown  In  black  and  white  on  the  display 
monitor. 

2,  Based  on  its  gray  level,  each  pixel  of  the 
Image  is  assigned  a  corresponding  pseudo¬ 
color  for  easy  viewing. 
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3.  A  Sobel  edge  operator  is  moved  sequentially 
through  the  entire  image  for  edge  enhance¬ 
ment.  Appendix  A  describes  the  Sobel  edge 
operator . 

4.  In  order  to  eliminate  the  noise  produced  by 
the  edge  operation  and  the  noise  appearing 
on  the  original  image,  a  lowpass  filter  is 
passed  through  the  whole  image.  The  low- 
pass  filter  is  described  in  Appendix  B. 


This  process  of  comparison  continues 
sequentially  for  all  unlabeled  pixels  in 
the  image.  When  the  process  is  completed 
up  to  this  point,  the  SUM^^  is  divided 
by  Cp  to  obtain  an  average  pixel  value 
for  a  particular  region  category  n  as 
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5.  After  the  lowpass  filtering,  a  technique 
called  "region  growing"  is  employed  to 
merge  together  pixels  which  have  similar 
gray  values.  A  commonly  used  simple  re¬ 
gion  growing  technique  Is  explained  in 
Appendix  C.  For  out  caee»  the  region  grow¬ 
ing  technique  was  modified  to  become  a  two- 
pass  process  as  described  below.  First  a 
threshold  value  T  is  sec.  The  selection 
of  a  threshold  value  T  requires  a  lengthy 
process  (Nagao  and  Matsuyama,  1980)  which 
will  be  described  in  detail  in  a  separate 
subsection  of  this  paper. 

Once  Che  threshold  value  T  is  selected, 
the  first-pass  of  the  region  growing 
operation  is  performed  as  follows.  A 
control  pixel  for  region  n  (or  category 
n) ,  called  selected  arbitrarily 

from  Che  image.  Usually  the  unlabeled 
upper  left  pixel  in  the  image  is  assigned 
as  the  control  pixel  "^^e  next  step 

is  CO  sequentially  compare  the  gray  value 
of  each  pixel  with  that  of  P..  as  given 
by  (1): 

where  G(l,j)  Is  the  gray  value  of  the 
pixel  P(l,j)  and  is  the  gray  value  ol 
Che  control  pixel  P^.^.  The  subscript  n 
signifies  chat  the  control  pixel  is  for 
Che  category  n. 

The  gray  values  of  all  pixels  that  meet 
the  inequality  (1)  are  Chen  summed  and  the 
result  la  added  to  the  gray  value  o(  the 
control  iMxvl.  VUv  ri'miUinj,  ,i,„„ 
saved  in  a  specially  designed  memory.  For 
a  particular  category  n  this  sununed  quan¬ 
tity  is  expressed  as  SL'M^,  and  is  given 
below: 


"  2  Cd.j).  i2) 

ail  pixels  met 
by  irequality  «L) 

Otherwise,  the  pixel  under  examination  is 
left  unlabeled.  At  the  same  time  a 
counter,  is  Incremented  by  1  when  a 

pixel  is  added  to  the  This  counter 

starts  with  a  content  of  i  so  that  the 
final  count  of  the  counter  will  Indicate 
the  total  number  of  pixels  added  to 
the  SUMn*  This  Is  expressed  for  the  cate¬ 
gory  n  as  follows: 


1  -►  Number  of  pixels  met 

by  Che  inequality  in  (1)  (3) 


The  next  step  is  to  repeat  the  entire  pro¬ 
cess  described  above,  or  computing  (1) 
through  (4)  for  the  remainder  of  the 
pixels  left  unlabeled  until  all  pixels  on 
Che  image  are  labeled*  end  eech  pixel 
belongs  to  a  particular  region.  This 
completes  the  first-pass  of  the  region 
growing  process.  At  this  point,  a  number 
of  average  pixel  values  for  potential 
regions  will  be  computed. 

The  second-pass  of  the  region  growing  pro¬ 
cess  Is  similar  to  Chat  of  the  first-pass 
except  chat  the  control  pixel  values, 

are  now  replaced  by  Che  correspond¬ 
ing  average  pixel  values  Ap.  Also*  when 
the  value  of  a  pixel  is  compared  to  a 
particular  average  pixel  value,  and  if 
the  absolute  value  of  the  difference  of 
the  pixels  is  less  than  Che  threshold 
value,  the  pixel  under  examination  will 
be  merged  to  that  particular  average 
pixel,  rather  chan  adding  it  Into  Che 
quantity  SUM,^.  In  ocher  words,  the  value 
of  the  pixel  under  examination  will  be 
sec  equal  to  that  particular  average 
pixel  value  chat  It  Is  being  compared 
with.  This  process  Is  explained  in 
Appendix  C.  The  second  pass  of  Che 
region  growing  process  will  continue 
until  each  pixel  in  the  Image  belongs 
to  a  particular  region  category.  The 
entire  region  growing  process  will  then 
be  complete. 

6.  The  number  of  region  categories  created 
by  Che  application  of  the  region  growing 
process  usually  sxcseds  four.  A  slmpls 
pixel  grouping  routine  Is  then  used  to 
lurttier  group  pixels  in  the  region  grown 
Image  Into  exactly  four  categories.  This 
routine  functions  as  follows:  The  gray 
value  of  each  pixel  In  Che  region  grown 
image  Is  sequentially  examined.  If  it  is 
larger  than  or  equal  to  100,  it  Is  set  to 
150.  If  it  Is  less  than  100,  and  larger 
than  or  equal  to  45,  it  is  set  to  65.  If 
it  is  less  than  45,  and  larger  than  or 
equal  to  8,  it  is  set  to  25,  otherwise 
it  Is  sec  Co  4. 

After  the  pixel  grouping  process,  the 
entire  Image  Is  segmented  into  exactly 
four  categories,  with  each  category 
assuming  a  different  pixel  gray  value. 

The  first  three  categories  having  gray 
values  of  4,  25,  and  65  represent  water, 
fields,  and  forests,  while  Che  last 
category  having  several  boundaries  or 
edges  with  a  gray  value  of  150  belongs 
to  a  built-up  area. 

5.  Having  done  the  image  segmentation  for  a 
given  image  as  described  above,  Che 
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extraction  or  finding  of  water  regions 
from  the  image  is  very  straight-forward. 
Set  ail  the  pixels  having  a  gray  value 
of  ^  to  0i  and  the  remainder  of  pixels 
to  255.  Now,  we  have  an  inverted  binary 
Image  of  water  regions. 

9.  In  order  to  eliminate  noise  within  the 
water  regions  which  now  appears  as  white 
spots  of  gray  value  255  on  the  inverted 
image,  a  square  box  containing  a  1  by  1 
pixel  is  made  to  scan  through  the  entire 
image  from  the  top  left  corner  to  the 
bottom  right  corner  in  a  usual  way.  Any 
255  pixel  (or  255  pixels)  which  can  be 
successfully  surrounded  by  this  box  repre¬ 
sents  a  separated  noise  pixel  (or  pixels), 
and  thus  its  gray  value  will  be  replaced 
to  0.  The  size  of  the  square  box  will  be 
sequentially  Increased  from  1  by  1  pixel 
to  16  by  16  pixels,  and  for  each  Increase 
the  same  scanning  and  replacing  process 
will  be  repeated. 

10.  The  inverted  Image  is  now  Inverted  again, 
so  all  pixels  within  the  water  regions 
will  be  represented  by  a  gray  value  of 
255.  The  noise  on  the  non-water  regions 
will  appear  as  white  spots  of  gray  value 
255.  The  process  discussed  in  Item  9 
will  be  repeated  for  eliminating  noise 
from  the  non-water  regions. 

11.  The  water  regions  which  are  represented 
by  a  gray  value  255  will  be  sec  to  180. 

We  selected  the  gray  value  180  to  repre¬ 
sent  Che  water  regions  because  it  shows 
up  as  a  pleasant  dark  blue  color  in  the 
pseudo-color  domain.  This  completes  the 
entire  water  finding  processes. 

12.  A  number  of  selected  SAK  images  contained 
shallow  water  bodies.  The  water  finder 
described  failed  to  extract  water  regions 
from  these  images  because  Che  pixels 
within  the  water  regions  appeared  Coo 
bright  (have  higher  gray  values  compared 
to  that  of  those  pixels  within  a  regular 
water  region).  A  second  water  finder 
with  a  relatively  higher  threshold  for 
Che  region  growing  process  was  designed 
to  compensate  this  problem. 

13.  Dlscrtmtnatton  of  non-water  features 
having  water  like  gray  values  is  accom¬ 
plished  by  using  a  third  water  finder 
equipped  with  an  additional  capability. 
Just  like  Che  previous  two  water  finders, 
the  water  regions  and  non-water  regions 
having  water  like  gray  values  are 
extracted.  An  algorithm  called  connected 
components  (Hevenor  and  Chen,  1990)  is 
applied  next.  The  purpose  of  this  con¬ 
nected  components  routine  is  to  provide 

a  unique  label  for  each  extracted  region 
from  Che  previous  operation.  The  detail 
of  the  connected  components  is  described 
in  Appendix  D.  Various  properties  of 
each  region  such  as  area,  centroid, 
elongation,  perimeter,  compactness,  etc 
can  be  computed  for  discriminating  non- 
water  regions  from  water  regions.  For 
our  case,  onlv  the  property  area  needs  to 
be  computed  since  the  largest  water  region 
always  has  the  largest  area  among  ail 
regions  all  SAR  Images  tested. 


The  calculation  of  area  for  each  region 
can  be  made  by  taking  a  summation  of  all 
pixels  within  the  region.  Finally,  the 
gray  value  of  each  pixel  within  the 
region  having  the  largest  area  is  set  to 
180,  and  Che  remainder  of  regions  are 
eliminated  by  setting  their  pixels  to  0. 

Selection  of  an  Optimum  Threshold  Value: 

The  selection  of  a  threshold  value  for  the 
pixel  based  region  growing  plays  a  crucial 
role  in  the  entire  process.  The  selection  of 
the  threshold  value  should  be  adaptively 
determined  by  Che  image  data  under  analysis. 
Using  an  improperly  predetermined  fixed 
threshold  value  for  region  grow^g  would 
lead  to  a  serious  mistake  and  would  end 
up  in  either  undergrowing  or  overgrowing  of 
the  image.  For  our  experimentation,  the 
original  SAR  Image  of  512  by  512  pixels  was 
first  edge  enhanced,  and  _sii^oched  with_a 
lowpass  filter.  These  two  steps  tog.k_plac£_  ±i 
the  first  step  discussed  in  Appendix  E.  The 
smoothed  image  was  divided  into  64  blocks  of 
subimages  with  each  block  consisting  of  64 
by  64  pixels.  The  threshold  determination 
method  discussed  in  (Nagao  and  Matsuyama, 

1980)  was  applied  to  each  block,  and  the 
minimum  threshold  value  found  for  the  sub¬ 
images  was  selected  as  Che  threshold  value 
for  performing  region  growing  for  the 
entire  image.  The  threshold  determination 
method  discussed  in  (Nagao  and  Matsuyama, 

1980)  is  provided  in  Appendix  E,  The 
optimum  threshold  value  was  12  for  both  sets 
of  SAR  Images  taken  over  the  areas  of  Eliza¬ 
beth  City  and  Ashboro,  North  Carolina.  With 
this  threshold  value  the  majority  of  tmagee 
tested  were  region  grown  to  yield  four  catc- 
guries  of  urea  terrain  features  as  described. 

RESULTS  AND  DISCUSSION 

The  software  for  Che  three  automated  water 
finders  was  successfully  tested  using  a 
group  of  18  images  selected  from  two  sets  of 
SAR  Imagery.  These  two  sets  of  SAK  imagery 
were  discussed  previously  in  the  section 
System  Description.  Test  results  for  Che 
water  finders  are  summarized  In  Table  1.  The 
name  of  the  specific  finder  (or  finders)  used 
In  each  Image  is  also  indicated  in  Table  1. 

For  illustration  purposes,  only  the  results 
obtained  from  three  Images  will  be  presented. 
It  is  seen  from  Table  I  that  for  most  test 
lisages  a  good  to  fair  result  was  obtained. 

Only  two  images  are  found  to  be  marginally 
acceptable.  The  following  observations  can 
be  made  from  the  results: 


Table  1.  Test  Results  for  Water  Finders 
Using  Two  Secs  of  SAR  Images 


Image 

Names 

Figure 

Number 

Finders 

Used 

Results 

UnfOO? 

NS 

Water 

Good 

Unf014 

NS 

Water 

Good 

Unf026 

la 

Water 

Good 

Unf033 

2a 

Waterl 

Fair 

Unf034 

NS 

WsCerl 

Fair 

Unf051 

NS 

Waterl 

Fair 

Unf052 

NS 

Waterl 

Marginal 

Unf053 

NS 

WaCer/WaCer2 

Fair 

l*nf056 

NS 

WaCer2 

Good 
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UnflXl 

3a 

Water2 

Good 

UnfllB 

NS 

Watex/Water2 

Good 

Unfl43 

NS 

Waterl 

Marginal 

Uii£163 

NS 

Waterl 

Good 

Unfl65 

NS 

Water2 

Fair 

Unfl67 

NS 

Waterl 

Fair 

Unfl82 

NS 

Water! 

Fair 

Erlml3 

NS 

Water2 

Fair 

Erlm27d 

NS 

Water2 

Good 

NS  •  Not  Shovm 


1.  Very  small  or  narrow  water  regions*  such 
as  the  narrow  stream  on  the  lower  left 
corner  of  the  first  original  image* 
Figure  la*  disappeared  after  the  edge 
operation  and  region  growing  processes. 
This  is  Illustrated  In  the  resultant 
Figure  lb.  The  first  finder,  named 
Water,  was  used  for  this  test. 


2.  To  extract ' shallow  water  regions  which 
have  relatively  high  pixel  gray  values* 
a  second  water  finder*  named  Waterl*  with 
a  higher  threshold  value  of  14  was  used. 
The  second  original  Image*  Figure  2d*  and 
the  corresponding  result*  Figure  2b, 
show  this  point.  The  narrow  and  lone 
triangle  appearing  on  Che  right  side 
of  Figure  2a  was  due  to  a  misalignment 
of  the  film  while  in  digitization.  The 
triangle  clearly  shows  an  edge  of  the 
film.  This  mistake  caused  an  error  In 
water  finding  which  appears  as  two 
near  vertical  and  two  horizontal  black 
lines  on  the  right  side  of  Figure  2b. 
Several  email  water  reglone  reeulclng 
on  Che  left  side  of  Figure  2b  are 
errors  due  to  an  Imperfect  region 
growing. 


Figure  le.  The  Original  SAR  Image*  Un£026.  Figure  2e.  The  Original  SAR  Image*  Uaf033, 


Figure  lb. 


The  Result  of  Applying  the  Water 
Finder  to  the  SAR  Image*  Unf026. 


Figure  2b. 


The  Result  of  Applying  the  Waterl 
Finder  to  the  SAR  Image .  Unf033. 


Discrimination  of  non-water  dark  regions 
whose  pixels  have  similar  gray  values  to 
those  in  water  regions  appears  to  be 
very  effective.  The  airfield  runways 
that  appeared  on  the  third  original 
image,  Figure  3a,  are  the  example  of 
this  discrimination.  They  are  completely 
removed  as  shown  in  the  resultant  Figure 
3b.  A  third  finder.  Wacer2.  was  used  to 
accomplish  this  discrimination. 


Figure  3a.  The  Original  SAR  Image*  Unflll. 


Figure  3b.  The  Result  of  Applying  the  Water2 
Finder  to  the  SAR  Image,  Unflll. 


CONCLUSIONS 


3.  Shallow  water  regions,  which  appear  con¬ 
siderably  brighter  than  ordinary  water 
regions  on  an  SAR  image*  can  be  extrac*' 
ted  by  using  a  separate  water  finder 
with  a  higher  threshold  value  for  the 
region  growing  process. 

4.  Discrimination  of  non-water  features 
which  have  approximately  the  same  gray 
values  as  that  of  water  regions  can  be 

a c c omp 11s hed  by^addlng  two  more  aigorTt hms 
(connected  components  and  region  pro¬ 
perty  computation)  on  the  first  water 
finder  as  described.  The  typical  non¬ 
water  terrain  features  of  this  nature 
Include  airfield  runways  and  shadows. 
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APPENDIX  A  SOBEL  EDGE  OPERATOR 


The  SobeL  operator  is  a  3  by  3  pixel  non¬ 
linear  edge  enhancement  mask  which  is  multi¬ 
plied  sequentially  with  all  pixel  values  in 
an  image  to  produce  a  pattern  of  more  pro¬ 
nounced  edges. 

The  weights  for  the  Sobel  mask  are  shown 
below: 

-10  1  12  1 
-2  0  2  0  0  0 

-10  1  -1  -2  -1 


Automated  finding  of  water  regions  from 
SAR  Imagery  can  be  effectively  accom¬ 
plished  by  properly  applying  a  set  of 
Image  processing  and  computer  vision 
algorithms  sequentially. 

Automated  water  finders  can  be  simply 
implemented  by  using  the  concept  of 
the  automated  terrain  segmenter*  which 
was  previously  developed,  as  a  pre¬ 
processor. 


x-direct ion  y-dlrectlon 

Assume  a  block  of  3  by  3  pixels  to  be  mulcl- 
plied  with  the  Sobel  mask  centered  at  the 
point  (l,j)  and  having  a  gray-value  dis¬ 
tribution  as  given  below: 
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APPENDIX  B  LOWPASS  FILTERS 

The  most  commonly  used  lowpass  filter  con¬ 
sists  of  a  mask  of  3  by  3  pixels.  It  is 
multipled  sequentially  with  all  pixel  values 
on  a  given  image  to  produce  a  smooched 
pattern. 


The  weights  for  the 
below: 

lowpass 

mask  are  shown 
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1 

1 

1 

1 

1 

1 
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Assume  a  block  of  3  by  3  pixels  to  be 
multiplied  with  Che  lowpass  mask  centered 
at  the  point  (l,j)  and  having  a  gray— value 
distribution  as  given  below: 

^0 

A. 

At 

A? 

Ffi.jj 

As; 

A  = 

Au 

Then,  the  magnitude  of  the  resultant 
which  will  replace  F(l,j).  will  be 

k-0 

For  more  sophisticated  lowpass  niters,  the 
size  of  the  mask  can  be  increased  to  5  by  5 
or  7  by  7  pixels,  and  have  values  consist¬ 
ing  of  all  I's.  However,  the  processing 
time  of  using  these  filters  will  be  Increased 
accordingly . 

APPENDIX  C  SIMPLE  REGION  GROWING 

The  simple  region  growing  method,  based  on 
pixel  gray  value,  consists  of  the  following 
steps : 

Step  1:  If  all  pixels  in  a  given  image  are 
labeled,  then  end.  Else  take  an 
unlabeied  pixel  and  assign  a  new 
unused  region  number. 

Step  2:  If  the  absolute  diff  erence  of  gray 
value  between  the  new  labeled  pixel 
and  Its  neighboring  pixels  are  less 
chan  the  threshold  value,  merge  the 
neighboring  pixels  and  assign  them 
the  same  region  number. 

Step  3;  Iterate  Step  2  until  no  pixels 

adjacent  to  the  new  labeled  region 
can  be  merged. 


Step  4:  Go  to  Step  1. 


appendix  d  connected  components 

The  purpose  of  the  connected  components 
routine  is  to  provide  a  unique  label  for  each 
component  of  1  pixel  in  the  binary  Image 
that  represents  water  regions.  Each  label 
is  a  number  chat  is  assigned  to  every  pixel 
In  a  given  connected  component.  This  label¬ 
ing  operation  can  be  performed  by  scanning 
Che  entire  binary  image  with  a  3-pixel  by 
3-pixel  array  and  considering  the  following 
pattern: 

C  B  E 

D  A 

If  we  scan  along  the  image  from  left  to 
right  and  from  Cop  Co  bottom,  and  if  pixel 
A  Is  the  pixel  presently  being  considered 
and  it  has  a  value  of  1,  then  a  label  must 
be  assigned  to  A.  The  pixel  at  D,  C,  B  and 
E  have  already  been  labeled  and  can  be  used 
to  label  A.  If  the  pixels  at  D,  C,  B  and 
£  are  all  0,  then  A  is  given  a  new  label.  If 
pixels  C,  B  and  £  are  all  0  and  D  ■  1,  Chen 
A  is  given  the  label  of  D.  Each  possible 
construction  of  O's  and  I's  for  the  pixels 
D,  C,  B  and  E  must  be  considered  when  pro¬ 
viding  a  label  for  A.  If  two  or  more  pixels 
In  Che  sec  D,  C,  B  and  £  are  equal  to  1  and 
they  all  have  the  same  label,  then  A  Is  also 
given  Che  same  label.  The  real  difficulty 
comes  when  two  or  more  of  the  pixels  D,  C, 

B  and  £  have  different  labels.  This  can 
occur  when  two  or  more  separate  components, 
which  were  originally  assigned  different 
labels  are  found  to  be  connected  at  or  near 
pixel  A.  For  these  cases  Che  pixel  A  is 
given  the  label  of  any  one  of  the  pixels  D, 

C,  B  or  E,  which  has  a  value  equal  to  1.  An 
equivalence  list  consisting  of  pairs  of 
equivalent  labels  Is  formed.  After  the  binary 
Image  has  been  completely  scanned,  the 
equivalauce  list  is  restructured  in  Ruch  a 
way  as  to  contain  a  number  of  lists.  Each 
of  these  lists  contains  all  of  the  equiva¬ 
lent  labels  associated  with  a  particular 
connected  component.  A  new  label  is  Chen 
given  Co  each  of  the  new  Hats  and  Is 
assigned  to  each  of  the  appropriate  pixels. 

APPENDIX  E  THRESHOLD  DETERMINATION  METHOD 

The  following  is  the  adaptive  threshold 
determination  algorithm  used  for  our 
experimentation. 

Step  I:  Differentiate  the  smooched  image 
using  the  operator 

d(l.J)  -  max  |G(i.J)  -  G(l+k,J+m) |  , 
-Kkil 
-limSl 


where  G(1,J)  and  denote  the 

gray  value  and  the  diff erentlal- value 
at  a  point  (1,  J) . 
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Step  2:  Divide  the  differentiated  image 
into  M  blocks  of  64  by  64  pixel 
subimages  and  make  a  histogram 
hj^(d)  of  the  differential  values 
d(i,j)  in  the  n-th  block  of  sub- 
image  (n  »  1,2 . M) . 

Step  3:  The  Valley-Detection  Algorithm. 

For  each  histogram  hp(,d),  find 
the  minimum  value,  d p  which 
satisfies  the  following  inequali¬ 
ties: 

h,,(d„)  <  k)  fork  =1.2 . N, 

where  d*  denotes  the  differential 
value  for  which  histogram 
has  the  maximum  population.  N  is 
set  initially  to  9.  The  value  of  N 
will  be  reduced  from  9  to  8,  and 
from  8  to  7  sequentially  until 
a  satisfactory  result  Is  obtained. 

Seep  4:  Find  the  minimum  value  among  the 
dp  for  all  blocks  of  subimages, 
and  make  it  the  threshold  value 
T  for  all  areas  of  the  image. 

That  is 


T  -  min  d 


